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Abstract 

We show that, in nonequilibrium systems with small heat flows, there is a 
time-scale dependent effective temperature which plays the same role as the 
thermodynamical temperature, in that it controls the direction of heat flows 
and acts as a criterion for thermalization. We simultaneously treat the case 
of stationary systems with weak stirring and of glassy systems that age after 
cooling and show that they exhibit very similar behavior provided that time 
dependences are expressed in terms of the correlations of the system. We 
substantiate our claims with examples taken from solvable models with non- 
trivial low-temperature dynamics, but argue that they have a much wider 
range of validity. We suggest experimental checks of these ideas. 

75.40.Gb, 75.10.Nr, 02.50.-r, 05.20.-y 
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I. INTRODUCTION 



No physical system is ever in thermodynamical equilibrium. When we apply thermody- 
namics or statistical mechanics, we idealize the situation by assuming that "fast" processes 
have taken place, and "slow" ones will not: hence, we define an observation time scale which 
distinguishes these two kinds of processes ||. It follows that the same system can be at 
equilibrium on one scale, and out of equilibrium on another, and, more strikingly, that it 
can be at equilibrium, but exhibiting different properties, on two scales at once. 

Since the assumption of thermal equilibrium lies at the heart of statistical mechanics, 
it is usually hard to make these considerations without a strong appeal to one's intuition. 
We show in the following that they can be made in fact quite precise for a class of systems 
characterized by very slow energy flows. These systems are out of equilibrium, either because 
they are very gently "stirred", i.e., work is constantly made on them, or because they have 
undergone a quench from higher temperatures a long time ago. 

The most typical example of such a system is a piece of glass that has been in a room 
at constant temperature for several months. Since the glass itself is not in equilibrium, we 
have in principle no right of talking about "the temperature of the glass", but only about 
the temperature of the room. However, we may legitimately ask what temperature would 
indicate a thermometer brought into contact with the glass and, again, we would be very 
surprised if it did not coincide with the room temperature. We would be even more surprised 
if, putting two points of the glass in contact with both ends of a copper wire, a heat flow 
were established through it. 

In other words, although we know that equilibrium thermodynamics does not apply 
for the glass, we implicitly assume that some concepts that apply for equilibrium are still 
relevant for it. This is not because the glass is "near equilibrium" but rather because it has 
been relaxing for a long time, and therefore thermal flows are small. 

Many attempts have been done to extend the concepts of thermodynamics to nonequi- 
librium systems — such as systems exhibiting spatiotemporal chaos or weak turbulence P,[7l . 
In this context, Hohenberg and Shraiman have defined an effective "temperature" for 
stationary nonequilibrium systems through an expression involving the response, the corre- 
lation and the temperature of the bath. A closely related expression appears naturally in 
the theory of nonstationary systems exhibiting aging such as glasses. 

We show here that this expression indeed deserves the name of temperature, because (i) 
the effective temperature associated with a time scale is the one measured on the system by 
a thermometer, in contact with the system, whose reaction time is equal to the time scale, 
(ii) it determines the direction of heat flows within a time scale and (iii) it acts as a criterion 
for thermalization. 

We shall here consider simultaneously two different conditions in which a regime with 
small flows of energy exists: 

1. Ordinary thermodynamical systems in contact with a heat bath at temperature T that 
are slowly driven ( "stirred" ) mechanically. The driving force is proportional to a small 
number which we shall denote D. 

The observation time and the manner of the stirring are such that, for long enough 
times, the system enters a stationary, time-translational invariant (TTI) regime [ID|: 
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one-time average quantities are independent of time, two-time quantities depend only 
upon time-differences, etc. Stationarity is a weaker condition than thermodynamical 
equilibrium, since it implies loss of memory of the initial condition but not all other 
properties that are linked with the Gibbs-Boltzmann distribution. 

2. Purely relaxational systems that have been prepared through some cooling procedure 
ending at time t — and are kept in contact with a heat bath at a constant temperature 
T up to a (long) waiting time t w (as in the example of the glass). t w is also usually 
called "annealing time" in the glass literature. In this case physical quantities need 
not be TTI, and in interesting cases they will keep a dependence upon t w (and also, 
in many cases, upon the cooling procedure) for all later times t — r + t w . 

We shall treat in parallel the "weak stirring" (D — > 0) and the "old age" (t w — > oo) 
limits: both taken after the thermodynamical limit of infinite number of degrees of freedom. 
We show that they lead to the same behavior, from the point of view of thermalization 
and effective temperatures, provided that one expresses time dependences in terms of the 
correlations of the system [[11]]. 

As a test of our ideas, we discuss thermalization in the context of the mean field theory of 
disordered systems, or the low-temperature generalization of the Mode Coupling Equations, 
but the nature of our results makes us confident that they have a much wider range of 
validity. 

In Section II we recall the generalization of the fluctuation-dissipation relation to the 
nonequilibrium case. In Section III we consider the reading of a thermometer coupled 
to a system: when the system is in equilibrium we show that the thermometer measures 
the temperature of the heat bath, while when it is out of equilibrium it measures different 
effective temperatures depending on the observation time scale. These effective temperatures 
are equal or higher than the one of the bath and are closely related to the FDT violation 
factor [1r^,jS|j9| introduced to describe the out-of-equilibrium dynamics of glassy systems. In 
Section IV we recall how time scales or correlation scales are defined in systems with slow 
dynamics. We then argue that, if the FDT violation factor is well defined within a time 
scale, a single degree of freedom thermalizes within that time scale to the corresponding 
effective temperature. In Section V we extend this analysis to several degrees of freedom 
and show that the effective temperature determines the direction of heat flows, and can be 
used as a thermalization criterion. In Section VI we discuss various phenomenological "Active 
temperature" ideas that have been used for a long time in the theory of structural glasses. 
Our conclusions are summarized in Section VII, where some experimental implications of 
our work are suggested. 

II. THE FLUCTUATION-DISSIPATION RELATION OUT OF EQUILIBRIUM 

Let us consider a system with N degrees of freedom (si, Sjy), whose dynamics is 
described by Langevin equations of the form 

Si = bi(s) +m(t) , (2.1) 
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where t]i(t) is the Gaussian thermal noise. For unstirred systems, we consider purely re- 
laxational dynamics, where the average velocity 6j(s) is proportional to the gradient of the 
Hamiltonian E(s): 

fc(*) = -£r<^. (2.2) 

The symmetric matrix V is related to the correlation function of the noise rj by the Einstein 
relations 

(r H {t)r h {t')) = 2TT ij 8{t-t') , (2.3) 

where T is the temperature of the heat bath. Averages over the thermal history, i.e., averages 
over many realizations of the same experiment with different realizations of the heat bath, 
will be denoted by angular brackets. We assume of course (r]i(t)) = 0, Vi,t. We have 
chosen the temperature units so that Boltzmann's constant is equal to 1. The equilibrium 
distribution is then proportional to the Boltzmann factor exp(—E/T). 

For stirred systems, we add to hi a perturbation proportional to D, that cannot be 
represented as the gradient of a function (i.e., is not purely relaxational), but is otherwise 
generic. We then have W = (J2ibi(s)si) > at stationarity, meaning that work is being 



made on the system [|Tg 



We denote the observables (energy, density, magnetization, etc.) by O(s). Throughout 
this work we shall denote by t w or t' the earliest time (to be related to the waiting time), t 
the latest time, and r the relative time t — t w . These times are measured, in the case of the 
unstirred systems which exhibit aging, from the end of the cooling procedure. 

Given two observables 0\ and 2 , we define their correlations Ci 2 (t,t w ) = 
(Oi(t)0 2 (t w )) — (Oi(t))(0 2 (t w )), and their mutual response 

#12 (Mw) = TTjQ > (2-4) 
Sh 2 {t w ) 

where h 2 appears in a perturbation of the Hamiltonian of the form E — > E — h 2 (t)0 2 . 
Obviously, causality implies Ri 2 (t, t w ) = for t < t w . It is also useful to introduce the 
integrated response (susceptibility) 

X 12 (Mw) = f dt'R 12 (t,t') . (2.5) 

Let us now make a parametric plot of x(£, ^w) vs - C(t, t w ) for several increasing values of 
t w . We thus obtain a limit curve lim^^oo x(t, t w ) = x(C). In the case of a weakly driven 
system, we wait for stationarity and plot x(t — t w ) vs. C(t — t w ) for several decreasing values 
of the driving D. We thus obtain the curve liniD^o x(t ~ ^w> D) = x(C). 

The fluctuation-dissipation theorem (FDT) relates the response and correlation function 
at equilibrium. One has 

Ri2{t -t w ) = - — , (2.6a) 

Xn(t - t w ) = ^ (C 12 (0) - C 12 (t - t w )) . (2.6b) 
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FIG. 1. The susceptibility x{t — t w ) vs. the auto-correlation function C(t — t w ) for the model of 
Appendix [A] once stationarity is achieved. The parameter D is equal to 0.05, 0.375, 0.025, 0.0125 
respectively, from bottom to top. The dots represent the analytical solution for the limit D — > 0. 
One sees that, in this limit, the FDT violation factor X{C) tends continuously to the dotted 
straight lines. The value of C at the breakpoint is C , the Edwards-Anderson order parameter 
or the ergodicity breaking parameter in the language of the MCT. 

If the equilibrium distribution is asymptotically reached for t w — ► oo (or D — ► in the case 
of strirred systems) the FDT implies that the limit curve x{C) is a straight line of slope 
-l/T. 

However, there is a family of systems for which the limiting curve x(C) does not approach 
a straight line (see Figs. and [|). For driven systems this means that the slightest stirring is 
sufficient to produce a large departure from equilibrium even at stationarity (Fig. [I]), while 
in the case of a relaxational system it means that the system is unable to equilibrate within 
experimental times (Fig. Q). 

Let us denote by —X(C)/T the slope of the curve x(C): 

d x (C) _ X(C) 



dC 



(2.7) 



This corresponds to 



*U = f (2.8) 



where the derivative is taken with respect to the earlier time. We have thus defined X(C), 
the FDT violation factor for nonequilibrium systems with slow dynamics. 

When FDT holds, for D — > (or t > t w — > oo respectively) we can treat the system as 
being in equilibrium and X tends to 1 in the limit. When this does not happen, we may 
inquire about the physical meaning of X(C). In order to answer this question, let us first 
recall the relationship there is between the FDT and the equipartition of energy. 
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FIG. 2. The susceptibility x(t,tw) vs. the auto-correlation function C(t,t w ) for the model of 
Appendix |A| at T < T g (D is strictly zero). The full curves correspond to different total times t, 
equal, from bottom to top, to 12.5, 25, 37.5, 50,75 respectively (i w > i/4 throughout). The dots 
represent the analytical solution when t w — » oo. Neither x(f, i w ) nor C(t,t w ) achieve stationarity 
(see Fig. |3| below). 

III. FREQUENCY DEPENDENT THERMOMETERS THAT MEASURE 

EFFECTIVE TEMPERATURES 

We use a harmonic oscillator of frequency u a to measure the "temperature" of a degree 
of freedom O(s) (O(s) may be the energy, or some spatial Fourier component of the mag- 
netization). At the waiting time t w we weakly couple the oscillator to the system via O(s), 
while we keep the system in contact with a heat bath at temperature T. We wait for a 
short time until the average energy of the oscillator has stabilized. If the system were in 
equilibrium, by the principle of equipartition of energy we would have (E osc ) = T. 

Assuming linear coupling, the hamiltonian reads: 

-Etotai = E(s) + E osc + E int , (3.1) 

where 

£osc = h 2 + l - u 2 Q x 2 , (3.2) 

-E int = —aO(s)x . (3.3) 

The equation of motion of the oscillator reads 

x = -u 2 x + aO{t) . (3.4) 

In the presence of the coupling, if ax(t) is sufficiently small (an assumption we have to 
verify a posteriori) , we can use linear response theory to calculate the action on O of the 
oscillator: 
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Oft) = OJt) +a I dt' Ro(t, t')x(t') , (3.5) 
Jo 

where Ob(t) is the fluctuating term and where the response function Ro is defined by 

Ro ^ t) = J^{¥)- (3 ' 6) 

We assume moreover that the average {0 D (t)) exists (we set it to zero by a suitable shift of 
x) and that the fluctuations of Ob (in the absence of coupling) are correlated as 

(O h (t)O h (t w )) = C Q (t,t w ) , (3.7) 

where Co(t,t w ) is a quantity of O(N). 
The equation for x then reads 



x 



-ulx + aO h (t) + a 2 / dt' R (t, t')x(t') . (3.8) 
Jo 



Thus the oscillator takes up energy from the fluctuations of O, and dissipates it through the 
response of the system. Equation (3^8) is linear and easy to solve in the limit of small a 2 by 
Fourier-Laplace transform. One thus obtains the following results, whose proof is sketched 
in Appendix [B]. 

Consider first the case of stirred systems at stationarity, in which both the correlation 
and the response are TTI. The average potential energy of the oscillator reaches the limit 

2"o <* ) = - 2 (E osc ) = -^j- , (3.9) 

after a time ~ t c given by 

2lu. 



a 2 X"(u ) 



(3.10) 



We have defined 



POO 

Im / dt R(t) e il 
Jo 

roo 

C(lu) = Re / dtC(t)e luJt . (3.11) 



If we now define the temperature T {oj) as follows: 

T (lu) = (E osc ), (3.12) 

we obtain 

This is precisely the temperature defined by Hohenberg and Shraiman for the case of 
weak turbulence (the spatial dependence is encoded in O(s)). 
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Let us now turn to the case of relaxational dynamics where TTI is violated. Here we 
have to take into account the time t w at which the measurement is performed, and consider 
uj^ 1 <C t c t w . Then, a similar calculation yields for the energy of the oscillator 

{X K = ^(^osc)t w = „ , , (3.14) 

where the average (x 2 ) tw is taken on a comparatively short time stretch after £ w . This 
definition of the frequency- and waiting-time-dependent correlation Co(oo Q , i w ) and out-of- 
phase susceptibility x'o( u o,t w ) closely follows the actual experimental procedure for their 
measure: one considers a time window around t w consisting of a few cycles (so that phase 
and amplitude can be defined) and small enough respect to t w so that the measure is "as 
local as possible in time". In fact, these two quantities are standard in the experimental 
investigation of aging phenomena in spin-glasses [p0| -[22[ . 



The natural definition of the frequency- and time-dependent temperature of O is then 

ToK, fw )= X ,, M ■ (3-15) 

If equilibrium is achieved, then the temperature is independent of £ w , of the frequency uj 
and of the observable O, and coincides with that of the heat bath. The index O recalls 
that the effective temperature may depend on the observable. The frequency- dependent 
temperature defined by either ( |3.13| ) or (|3.15|) is compatible with the Fourier transformed 



expression of the FDT violation factor (|2.8f) provided that it does not vary too fast with u . 
We show later that this is indeed the case for systems with slow dynamics. 

For this definition we have chosen somehow arbitrarily an oscillator as our thermometer. 
However, we show in Appendix |CJ that the same role can be played by any small but 
macroscopic thermometer, weakly coupled to the system. The role of the characteristic 
frequency uj is then played by the inverse of the typical response time of the thermometer. 

Now To(ui,t w ) only deserves to be called a "temperature" if it controls the direction of 
heat flow. A first way to check whether this is the case is to consider an experiment in 
which we connect the oscillator to an observable 0\ and let it equilibrate at the temperature 
T 0l (w,t w ). We then disconnect it and connect it to another observable 2 and let it equi- 
librate at the temperature To 2 {uj, t w ). This is not like a Maxwell demon, since the times of 
connection and disconnection are unrelated to the microscopic behavior of the system. The 
net result is that an amount of energy T 0i {oj, t w ) — T . 2 {oj, t w ) has been transferred from the 
degrees of freedom associated with 0\ to those associated with 2 : therefore the flow goes 
from high to low temperatures. This is an actual realization of the idea of "touching two 
points of the glass with a copper wire" described in the Introduction. 

This observation also suggests a possible explanation of the fact that all FDT violation 
factors which we know of are smaller than one: if Tq{uj, t w ) were smaller than the tempera- 
ture of the heat bath, it could be possible, in principle, to extract energy from the bath by 
connecting between it and our system a small Carnot engine. The argument can be made 
sharper by considering a stationary situation in a weakly stirred system: but to argue that 
this situation would lead to a violation of the Second Principle one needs to prove that the 
equilibration times of the Carnot engine are short enough to make the power it produces 
larger than the dissipated one. 
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IV. EQUILIBRATION WITHIN A TIME-SCALE 



Before discussing the effective temperature as an equilibration factor we need to introduce 
some general features of the time-evolution of systems with slow dynamics. We first define 
the time-correlation scales and we then argue that, if the FDT violation factor is well 
defined within a time scale, a single degree of freedom thermalizes within that time scale at 
the corresponding effective temperature. 

A. Correlation and response scales 

Systems having a long-time out-of-equilibrium dynamics tend to have different be- 
haviours in different time-scales. Let us start by describing them for long-time dynamics in 
the purely relaxational case. In Appendix we give a formal definition of correlation scale, 
following Ref. §. 

Consider first, as an example, the dynamics of domain growth |23| for the Ising model 
at low but non-zero temperatures. The autocorrelation function C(t,t w ) = jrJ2i Si(t)si(t w ) 
exhibits two regimes: 

• At long times t, t w , such that t — t w -C t w , the correlation function shows a fast decay 
from 1 (at equal times) to m 2 , where m is the magnetization. This regime describes 
the fast relaxation of the spins within the bulk of each domain. 

• At long and well separated times (t — i w ~ t w ), the correlation behaves as a function 
of L(t w ) I L(t) , where L{t) is some measure of the typical domain size at time t. 

We refer to these two scales as "quasiequilibrium" and "coarsening" (or "aging"), respec- 
tively. After a given time, the correlation rapidly decays to a plateau value m 2 , and the 
speed with which it falls below that value becomes smaller and smaller as t w grows. 

This example helps to stress the fact that different scales are defined as a function of 
both times (in this case t — t w finite, L(t w ) / L(t) finite) and are well separated only in the 
limit where both times are large. 

In Fourier space, this separation of scales is achieved by considering several frequencies 
uj and increasing times t w . Then, if we consider uj ~ const, for increasing t w , we probe the 
quasiequilibrium scale, while if we want to probe the aging or coarsening scale we have to 
consider smaller and smaller uj, keeping u)L(t w )/ L'(t w ) ~ constant. 

If we now consider two frequencies uj\ and uj 2 , and keep 

— ~ const. , (4.1) 

we shall probe, as t w — > oo, the same scale: it will be the quasiequilibrium scale if they both 
remain finite, or the coarsening (aging) one if we keep ui t 2L(t w ) / L' (t w ) ~ const. 

These considerations can be generalized to other systems with slow dynamics. In general 
there may be more than two relevant scales [p^,p|,p|,p4| , for example oj = const., ojtl/ 2 = 
const., ut w = const., etc. In any case, the condition for looking into the same scale via two 
successions of frequencies remains Eq. ([4.1|). 
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In Fig. ^ we show the numerical solution of our test model, defined in Appendix [A]. The 
auto-correlation function C(t,t w ) is plotted vs. the time difference t — t w for several waiting 
times in log-log scale. These plots (which are standard in the MonteCarlo simulations of 
spin-glasses f25|), show (i) that the system is out of equilibrium, since we have an explicit 
dependence on t w , and (ii) that there are at least two time-correlation scales. For short time 
differences t — t w the decay is fast, the auto-correlation function is TTI and it falls from 
1 to C EA . It then decays further from C EA to zero, more and more slowly as the waiting 
time increases. The quantity C EA is known in the language of spin-glasses as the Edwards- 
Anderson order parameter and in Mode Coupling Theory (MCT) as the non-ergodicity 
parameter: it measures the strength of the fast correlations {C EA = m 2 in domain growth). 
Quasiequilibrium and aging regimes are experimentally observed in real spin-glasses and 



polymer glasses [p6|,p2 



Let us now turn to a similar analysis for driven systems in the limit of weak driving 
energy D — > 0. In that case, even if the system reaches a TTI regime provided we wait long 
enough (the amount of waiting increases when the stirring rate D decreases), it sometimes 
happens that some correlations and responses acquire a non-trivial low frequency behaviour 



in the limit D — > P? |. For example, in Fig. |] we show the same plot as in Fig. |3], i.e., 
C(t — t w ) vs. t — t w , for different, but small, stirring rates. We see that there are at least 
two time scales: one, for short time differences, where the correlation decays rapidly from 1 
to C EA and one, for long time differences, in which it slowly decays from C EA to zero. The 
smaller the stirring, the slower the decay of C from C EA to zero. 

It is useful to consider, in these cases, frequencies that go to zero as some function of 
D. The condition that two sequences of frequencies uo 2 correspond to the same scale now 
reads 

lim — = const. . (4.2) 



B. Thermalization criterion for a single degree of freedom within a scale 

Let us now consider a system at a given waiting time t w or stirring rate D. Heat flows 
tend to zero as £ w — > 00 or D — > respectively. It is the reasonable to assume that, in these 
limits, the effective temperatures associated with any given observable O on a given time 
scale tend to equalize, provided that their limit is finite. We thus have, for example in the 
case of relaxational dynamics, for Ui/u 2 = const., 

lim T (a>i,£ w ) = lim T (u} 2 ,t w ) . (4.3) 

t w — >oo r w — >oo 

Similarly, we expect that each fixed frequency sooner or later thermalizes with the bath: 

lim T (u,t w )=T. (4.4) 
u)= const. 

In some important cases the effective temperatures defined by eqs. ( |3.13| ),( |3TT5| ) diverge in 
the limit t w — > 00 (D — > 0). In these cases the effective temperature T (uj, t w ) should diverge 
for the whole time scale. We shall dwell on this problem in Sec. C. 
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FIG. 3. The correlation function C(t,t w ) vs. t — t w . for the same model as in Fig. ||. From 
bottom to top t w = 20, 50, 100, 150, 200. The correlation decays rapidly from 1 to C EA ~ 0.73 
and then more slowly from C EA to 0. This second decay becomes slower and slower as t w increases. 
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FIG. 4. The correlation function C(t — t w ) vs. t — t w for the same model as in Fig. |T[ From 
bottom to top, D = 0.1, 0.075, 0.05, 0.025. The correlation decays rapidly from 1 to C EA ~ 0.73 
and then more slowly from C EA to 0. This second decay becomes slower and slower as D decreases. 
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Equations ( f4.3|) and ( |4.4j ) are trivially true in a system that reaches thermal equilibrium, 
where all effective temperatures eventually reach the temperature of the reservoir. How- 
ever, they also describe the situation in which smaller frequencies take longer to reach the 
temperature of the heat bath, in such a way that at each given (long) waiting time there 
are low enough frequencies that have a temperature substantially different (in all cases we 
know, higher) from that of the bath. In particular, Eq. (fOJ) allows us to answer a question 
we asked at the beginning: if a piece of glass has been kept at room temperature for several 
months, a thermometer whose response time is of order of a few seconds would measure 
the room temperature, but it would read a higher temperature if its response time is of the 
order of weeks. 

In fact, with the appropriate handling of time-scales || (see Appendix |D|), Eqs. ( |4.3|) 
and ( |4.4| ) make it possible to calculate the out-of-equilibrium relaxation of mean-field spin- 
glasses, and also the low temperature generalization of the mode-coupling equations for one 
single mode. We thus obtain a solvable example where Eqs. ( [4.3| ) and (O) hold. 



For the case of stirred systems, the D — > limit plays the same role as the t w — * oo limit 
in relaxational systems, provided that the time scales are suitably redefined as in Eq. ( (4.2| ). 
In particular, at each fixed value of uj one has lim^^o To(uj, D) = T, and, for U\/&2 = const, 
one has 

hm 2b (wi, D) = hm T (cj 2 , D) , (4.5) 

provided that the limits are finite. 

If one looks into a time scale within which the temperature is almost constant, one can 
indifferently use a temperature defined in terms of a frequency or a time. In the relaxational 
case one has 

To{u,t w ) =T {t,t w ) , (4.6) 
provided that t — t w ~ uj~ 1 . In the same situation one has 

^ w) _ 1 dC(t,t w ) (47) 
To(t,t w ) dt w 

Similar relations hold for the driven case. If, on the contrary, one considers values of fre- 
quency or time for which the temperature defined by Eq. fl4.7f) is not constant, one cannot 
directly relate it with the reading of a thermometer coupled to the system. 



C. Coarsening and the case of infinite effective temperature 

An important physical situation which deserves a special discussion is that of domain 
growth. In this case the effective temperature (associated, say, with the total magnetization) 
in the "coarsening scale", uL(t w ) /L'(t w ) = const., tends to infinity as t w — > oo. The curve x 
(magnetic susceptibility) vs. C (magnetization correlation) looks like in Fig. [|. This figure 
shows the results for a model that is equivalent to the O(N) ferromagnetic coarsening in 
three dimensions. The integrated response x becomes flat, as t w — > oo, in the aging regime. 
In other words, the long-term memory tends to disappear. 
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FIG. 5. The susceptibility x(Mw) vs. the auto-correlation function C(t, i w ) for the purely 
relaxational p = 2 spherical model (equivalent to the 0(N), N — > oo ferromagnetic coarsening in 
d = 3). The dots represent the analytical solution when i w — > oo. The total time i is equal to 20, 
50, 100, 200 (4* > t/4 throughout). 



Because experimental measures of aging are in general related to the response, this kind 
of system is sometimes referred to as exhibiting aging "in the correlations" (correlations are 
not TTI) but not in the response. One should stress, however, that the divergence of the 
effective temperature in the aging time scale can be extremely slow: for example, in the 
Fisher- Huse model for spin-glasses T (u,t w ) grows like a power of logt w . 

The fact that the effective temperature in the aging regime tends to infinity means that 
if we measure it at u small and fixed, the temperature T(u,t w ) grows with t w (while we 
are still probing the aging regime) and then starts falling, finally reaching the temperature 
of the bath. If we repeat the experiment with a smaller u, the overall behaviour is the 
same, but the highest temperature reached will be higher; and so on without an upper limit. 
The situation is similar, but experimentally more subtle, for "stirred" systems for which the 
effective temperatures diverge as D — > 0. In such cases we measure a system that has a 
given (small) driving D. If the system is such that for D — > the effective temperature 
of all but the fastest scale tend to infinity, we may never realize that this is the case in an 
experimental situation, as long as we are not able to repeat the experiment with smaller 
and smaller stirring. 

We show one such case in Fig. [| which represents the same coarsening problem as before, 
but made stationary by a "stirring" term in the equation of motion proportional to D. If 
one had performed an experiment on such a system for a fixed D, one would have observed 
high effective temperatures for separated times, and no evidence of thermalization. Only 
by letting D take smaller and smaller values one can notice that in fact these temperatures 
tend to infinity. 
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FIG. 6. The susceptibility x(t ~ *w) y s- the auto-correlation function C(t — t w ) for the p = 2 
spherical asymmetric model ("stirred" O(N), N — > oo ferromagnetic coarsening in d = 3) for 
different levels of asymmetry (stirring): from bottom to top, D = 0.8, 0.6, 0.4, 0.2. The dots 
represent the analytical solution in the limit of zero asymmetry. 



D. Systems with finite effective temperatures 

In mean-field models one can close the Schwinger-Dyson equations into a set of dynamical 
equations involving only the correlation and response functions. This allows us to solve them 
analytically for large times, and also to obtain their full solution, via a numerical integration, 
starting from a pair of initial times £ w = t — 0. We thus obtain the curve x{C) by integrating 
the response function and plotting it vs. C. Two cases in which there are only two time- 
correlation scales are the test model || considered here (cf. Appendix and the case of 
a particle moving in an infinite-dimensional random medium with short-range correlations 
p9| . Both the analytical and the numerical results exhibit thermalization within the aging 



regime. 

In the Sherrington-Kirkpatrick spin-glass model there is a full hierarchy of time scales 
and effective temperatures, a fact also confirmed numerically 0] and by Monte-Carlo sim- 
ulations ||. Another model with infinite many scales and with a full hierarchy of effective 
temperatures is that of a particle moving in an infinite-dimensional random medium with 
long-range correlations [51,53]. 



In Ref. pi |, the Monte-Carlo simulation of the "realistic" 3D Edwards- Anderson model 
for a spin-glasses was used to obtain the x{C) curves, which seem to approach a non-trivial 
curve for increasing t w . This suggests the existence of a hierarchy of time-correlation scales. 

The "trap model" for spin-glass dynamics [|32| violates Eq.(fOj) when an unusual choice 



of a parameter (a) is made. However, it is difficult to interpret the model, with this choice of 
parameter, as a phenomenological model stemming from a reasonable microscopic dynamics. 
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V. THERM ALIZATION OF DIFFERENT DEGREES OF FREEDOM WITHIN A 

TIME SCALE 



As we remarked in Section III, a bona fide temperature should control heat flow and 
thermalization. It is the aim of this Section to show that this is indeed the case for the 
effective temperature we have defined. If this were not the case, it would be possible to 
use our small oscillator to transfer heat from some degrees of freedom to others: in other 
words, by decorating our copper wire with a suitable frequency filter, we would observe heat 
flowing through it when it touches the two end of the glass. 

We shall argue that if different degrees of freedom effectively interact on a given time 
scale, then they thermalize on that scale. A useful — though not universal — criterion for 
"effective interaction" is that their mutual response (the response of one of the degrees 
of freedom to an oscillating field of the given frequency conjugate to the other degree of 
freedom) is of the same order as their self-response on that time scale. 

The argument is essentially the same both for relaxational systems that have evolved for 
a long time, or in stationary driven systems in the limit of small driving energy. We shall 
thus focus on the relaxational case only. 

We emphasize again that the thermalization of different degrees of freedom is well defined 
only in the limit of vanishing heat flow, i.e., long waiting times or vanishing stirring rates 
respectively. This is witnessed by the appearance of one (or more) well separated plateaus 
in the decay of the two-time correlation functions. 

Let us consider n modes (labelled by a, b = 1, . . . , n). One can write, in general, some 
Schwinger-Dyson equations for their correlations C = (C ab ) and responses R = (R ab )'- 



dC ab (t, t v 



9R a b(t, t 



— E ^ac{t) C c b{t, t w ) + 2T R ab (t w , t) + 

c 

+ £ / dt" D ac (t, t") R cb (t w , t") + ]T / dt" E ac (t, t") C cb (t", t w ) , (5.1a) 

c J c 

E V*c(t) Rct(t, t w ) + 5{t - t w ) 5 ab + E T dt " t") R cb(t", t w ) • (5.1b) 



As they stand, Eqs. ( |5.1| ) are just a way to hide our difficulties under the E, D carpet. 
Several approximations extensively used in the literature amount to various approximations 
of the kernels S and D. 

In equilibrium, the symmetries of the original problem allow us to write: 



Z ab (t-t w ) = ±^(t-t w ), (5.2a) 
T dt w 

R ab {t-U) = ^^{t-t w ). (5.2b) 

If we now make for S and D the approximation that are ordinary fuctions (instead of 
general functionals) of the correlations and the responses, we obtain the Mode- Coupling 
Approximation (MCA) 

D ab (C) = F ab (Q , S afc (C) = E F abjCd (QR cd , (5.3) 

c.d 



15 



where there is a model-dependent function F(q) such that 



OF d 2 F 
F ab (q) = -7T- ; F abM {q) = ~— ~— • (5-4) 

0<lab 0%bO<icd 

If the system equilibrates, one recovers the usual form of the Mode Coupling Theory 
(MCT) that is applied for supercooled liquids pifl . It will be instructive to recall how 
equilibrium is reached within the framework of eq. (|5.4|) . One recalls 

• The FDT conditions (^|); 

• Time translational invariance: functions of one time are just constant and functions 
of two times depend upon time differences only. 

• Reciprocity: C ab (t - t w ) = C ba (t - t w ). 

Putting in this information in Eqs. (|5.1|) and (|5.4|) one obtains 



dCab(t dt tw) = ~ E Ha* C cb (t - t w ) + 1 £ [D ac (0) C cb (t w - t w ) - D%C%\ 

+f E £ dt" D ac {t - t") dCcbi g t ~ tw) , (5.5) 

where D^.,C°^ stand for limits of widely separated times, i.e., 

C™ = lim C ac (t) , DZ = Van D ac {t) . (5.6) 

t — >oo t — ^oo 

This is the usual equilibrium MCT equation |34j . The equation for the response (|5.1b| ) with 
the same equilibrium assumptions becomes the time-derivative (divided by T) of Eq. ( |5.1a| ), 
as it should. 

In a case in which the system is unable to reach thermal equilibrium, like the low tem- 
perature phase of spin glasses, the solution of Eqs. (|5.1| ) will exhibit several time scales. 
The question as to how many scales one has to consider in order to close the dynamical 
equations can be answered for each model unambiguously by the construction in Ref. ||, 
suitably generalized to several modes. 

Here, for simplicity, we assume that there are only two relevant time scales (as in the 
coarsening example of the last section). We discuss later a model system that acts as an 
explicit example of this situation. We propose an ansatz for the long-time asymptotics, and 
then verify that it closes the equations ||. We assume (and later verify) that all two-time 
functions can be separated in two regimes: 

• Finite time differences with respect to the long waiting time, i.e., t — t w finite and 
positive, and t w — > oo. In this regime TTI and FDT hold. 

• Aging regime, corresponding to long and widely separated times, i.e., t ~ t w — > oo. In 
this regime neither TTI nor FDT hold. 
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For finite time differences and a long waiting time t w we have: 

Ca b DT (t - t w ) = lim C ab (t, t w ) , (5.7a) 



£ w — >oo 

vFDT, 



iCr (* - *w) = lim R ab (t, t w ) , (5.7b) 
with 

i XfiFDT 

Cf^EE lim W hm C™ T (t-t w ). (5.9) 

The aging regime is defined as the time-domain in which the correlations fall (more and 
more slowly) below C^ b A . For these times we denote 

C ab {t,t w ) = C ab (t,t w ) , R ab (t,t w ) = R a b(t,t w ) . (5.10) 

The separation ( |5.7| ) induces within the MCA a similar separation for E and D (cfr. 
Eqs. (|5-3 ) and (|5.4Q ), namely, for close times, 



D™ T (t - t w ) = lim D ab (t,t v 

c w — ^oo 



where FDT holds: 



Again, we can define 



Kb (t - U) = lim E ab (t, t w ) , (5.11) 

£ w — >oo 



i f)r>FDT 

Kb DT {t-t w ) = Y J ^{t-K)- (5-12) 



lim lim D F ab DT (t - t w ) ee Dg A , (5.13) 

£-tw — >oo r w — >oo 



and mark with the tilde the aging part of the kernels: 

D a b(t,t w ) = D ab (t,t w ) , S ab (t, t w ) = S a fe(t, t w ) . (5-14) 

In order to close the dynamical equations, we make an ansatz for the aging parts C ab (t, t w ) 
and R ab (t,t w ). For a problem in which the correlations vary only within two time scales, 
the natural generalization of the solution in Ref. f| is 

C a b(t,tw) = C ab (h ab (t w )/h ab (t)) , (5.15a) 
iWMw) = ^^(Mw) , (5.15b) 

where the X ab are constants and the h ab (t) are functions to be determined from the dynamical 
equations. Note that the derivative in ( |5.15|) is taken with respect to the earliest time. 

Remarkably, it turns out that one can close the equations with two different types of 
ansatz for the long-time aging behaviour. In Appendix [E| we show how this is done. In 
terms of the effective temperatures the meaning of these two possibilities is: 



17 



1. Thermalized aging regime. The effective temperatures associated with the observables 
0\, 2 are equal to each other for frequencies and waiting times in the aging regime: 
they are not necessarily equal to the temperature of the bath. At higher frequencies, 
they both coincide with the one of the bath. 

T\(uj, t w ) = T 2 (u, t w ) = T , t w -> 00 , C ah > C^ A , quasieq. . . 

Ti(u, t w ) = T 2 (u,t w ) , u -> , t w -> 00 , C ah < C^ A , aging. ' 

Not surprisingly, in this case we find that 0\ and 2 are strongly coupled (also in the 
aging regime) in the sense that the mutual responses 

Ru{t,tw) = T^U > R2i(t,t w ) = - , (5.17) 

where X > 0, are of the the same order of the self response functions. 

2. Unthermalized aging regime. The effective temperatures associated with the observ- 
ables 0\, 2 for combinations of frequencies and waiting times corresponding to the 
aging regime are neither equal to each other nor to that of the bath, while for higher 
frequencies they both coincide with the one of the bath. 

Ti(uj, t w ) = T 2 (u, t w ) = T , t w -> 00 , C ah > C^ A , quasieq. . . 

T\(uj,t w ) ^ T 2 (uj,t w ) ^ T , c^^O, t w -> 00 , C ab <C^ b A , aging. " ' 

In this case, 0\ and 2 are effectively uncoupled (in the aging regime), in the sense 
that: 

fli2(Mw)= T ^ » R 21 (t,t w )= T ^ , (5.19) 

where 

Xi 2 ,X 21 ->0. (5.20) 



We have not found any other way of closing the equations [|36|. We shall show below, 
in a particular case which can be numerically solved, that indeed either case (1) or case (2) 
take place. Let us remark that a solution with X\ 2 = X 2 \ 7^ 0, X\\ = X 22 = 1^0 and 
X 7^ Xx 2 is not compatble with the interpretation of X as an inverse temperature. We do 
not find such a solution in our test models and believe that it is not realizable in general. 

The considerations of this section can also be made in the limit of small stirring. 
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FIG. 7. The susceptibility Xn(Mw) and X22(t,t w ) vs. the corresponding auto-correlation func- 
tions Cn(t, t w ) and C22(t,t w ) for the uncoupled, aging systems. The slopes of the curves, i.e., the 
FDT violation factors and hence the effective temperatures are different. This corresponds to the 
unthermalized case. 
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FIG. 8. The susceptibilities Xij(t,t w ) vs. the correlation functions Cij(t,t w ) for the two aging 
systems of Fig. ^, this time weakly coupled. The FDT violation factors Xy(Cy) are almost parallel: 
thermalization is almost complete. 
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FIG. 9. The susceptibilities Xij{t ~ *w) y s- the correlation functions Cij(t — t w ) for two slightly 
driven system (D = 0.1), weakly coupled. The FDT violation factors Xij(Cij) given by the slopes 
of these curves are now the same: after a short transient corresponding to finite times all the curves 
become parallel. The aging regimes have thermalized. 

A. An explicit model with thermalization 

In order to test in a particular example that these asymptotic solutions are not only 
consistent, but are in fact reached, we solve numerically the mode coupling equations with 
two coupled modes, with F given by: 

F(q) = <f n + K 2 q p 22 . (5.21) 

We impose normalization at equal times of the autocorrelation of both modes: 

C 11 (t,t)=C 22 (t,t) = 1 . (5.22) 

As usual with the mode coupling equations p5j^8| - ^0| , ^| , one can find a disordered mean 
field model whose dynamics is exactly given by these equations. This is the model defined 
in Appendix |A|. 

Figures [7] and ^] show the numerical solution of the exact system of coupled integro- 
differential equations whose large-time asymptotic can be analytically obtained as in Ap- 
pendix [E|. In Fig. [7| and || we plot x(6) f° r two cases. In Fig. [7| we consider two 
uncoupled systems that evolve from the initial condition Cn(0, 0) = 6*22(0, 0) = 1 and 
Ci2(0,0) = 621(0,0) = 0. We see that the systems evolve independently and X\\ 7^ X22 
while X12 = X21 = (unthermalized case). In Fig. § we consider the same global system but 
now including a weak coupling between the two individual systems. Clearly, after a short 
transient associated to short times, all curves Xiji^ij) parallel. The systems aging-regime 
temperatures for the two subsystems have become equal (thermalized case). 
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It is interesting to note that, in this example, Onsager's reciprocity relations hold: R\ 2 = 
R21, C12 = C21. However, one can imagine situations in which they do not hold separately, 
but where however X 12 = X 2 \. 

For comparison, we show in Fig. ^|the corresponding plot for two weakly coupled systems 
in the limit of small stirring. 



B. Many scales 

Let us briefly discuss what would happen in the presence of many scales, each one with 
its own temperature. 

Consider, for instance, two observables 0\ and 2 , their associated autocorrelation func- 
tions Cn(t,t w ), C 22 (t, t w ), their integrated self- responses Xn(Mw), X22(t,t w ) and the ef- 
fective temperatures Ti(Cn), T 2 {C 22 ). We can plot Ti(Cu) vs. Cxi and T 2 (C 22 ) vs. C 22 . 
These two plots need not be the same, even if all scales are thermalized. Consider now a 
parametric plot of Cn(t,t w ) vs. C 22 (t, £ w ), in the limit of very large t w : it defines a function 
C11 = H(C 22 ) that allows one to calculate Cn for large t,t w , given C 22 (t, t w ). 

The condition for thermalization in every time-scale is then that the curve T 2 (C 22 ) coin- 
cides with the curve T 1 (7Y(C 22 )), both considered as functions of C 22 : 

T l {H{C 22 )) = T 2 {C 22 ) . (5.23) 

The deviation of T 1 (7i(C 22 )) from T 2 {C 22 ) is a measure of the degree in which the two 
observables are not thermalised. 

Actually, (|5.23|) was obtained as an ansatz for the dynamics of a manifold in a random 
medium within the Hartree approximation [JO], where the role of different observables is 
played by the displacements at different spatial wavelengths: 

Ck(t,t) = Hk(Ck=o) , 
T k (H k (C k=0 )) = T k=0 (C k=0 ) . (5.24) 

In this case, the modes at different k thermalize at the same effective temperature, although 
their mutual responses vanish, as implied by translational symmetry in mean. 



VI. COMPARISON WITH OTHER EFFECTIVE TEMPERATURES 



The idea of "Active temperatures' : 
developed since p3|-^5[ 



Tf in glasses goes back to the '40s and it has 
Here we recall it briefly, for the sake of comparison with the 
effective temperature that we have discussed. 

When cooling a liquid the time needed to establish equilibrium grows and, eventually, 
the structural change cannot keep pace with the rate of cooling: the system falls out of equi- 
librium and enters the glass transition region. It is then said that "the structure is frozen" at 
a temperature characterized by a Active temperature Tf. The Active temperature defined in 
terms of different quantities of interest, e.g., the enthalpy, the thermal expansion coefficient, 
etc., do not necessarily coincide. Furthermore, it has been experimentally observed that 
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glasses with the same Active temperatures arrived at through different preparation paths 
may have different molecular structures. The Active temperature is hence a phenomenolog- 
ical convenience and should not be associated with a definite molecular structure . 

The Active temperature is a function of the temperature of the bath T. At high temper- 
ature, when the sample is in the liquid phase, Tf = T. When the liquid enters the transition 
range Tf departs from T and Tf > T; and, finally, deep below the transition range, where 
the relaxation is fully stopped, Tf — > T g . The detailed bath-temperature dependence of the 
Active temperature in the region of interest is usually expressed by |f44| , |45| 



tf= -7(k^)" ((U) 

where the characteristic time r(Tf , T) depends both upon Tf (t) and upon the thermal history 
of the sample given by T(t). At equilibrium Tf = T and Tf = 0. This nonlinear difierential 
equation determines Tf(t) once one chooses r(Tf,T). A commonly used expression is the 
Narayanaswamy-Moynihan equation 



T (Tf, T) = r exp 



xA (l-x)A 



(6.2) 



where r G , A and x G [0, 1] are some constants. All the information about the dynamics of 
the system enters into Tf through these constants. 

In order to obtain the time relaxation of the quantity of interest, the picture is completed 
by proposing a given relaxation function, like the stretched exponential, the Davidson- 
Cole function, etc., and by introducing the Active temperature through the characteristic 
time r(T f ,T). See Ref. j|5| for an extensive discussion about the applications of T f to the 



description of experimental data. 

One may wonder whether is a relation between this fictive temperature of glass phe- 
nomenology and the effective temperature we have been discussing in this paper. 

First of all, one notes that the Active temperatures are deAned through the relaxation 
of the observables, unlike the one we consider here which are deAned rather in terms of 
Auctuations and responses. Both temperatures may depend upon the observable and upon 
the thermal history of the sample. 

The dependence of the efiective temperature deAned in Eq. (|3|) upon T depends on the 
model. In all cases FDT holds in the high temperature phase and To(cu, t w ) = T, the efiective 
temperature is equal to the bath temperature. When entering the low temperature phase, 
the temperature dependence of the effective temperature observed at fixed low frequency 
depends on the model. In certain simple mean-field (or low-temperature mode-coupling) 
models whose dynamics following a quench into the glass phase can be solved, one can 
compute To(u;,t w ) explicitly. Three different behaviours are found: 

• In the simple model we have been using as a test example in the previous sections ||, 
To(w,t w ) starts from To(uj,t w ) = T g at T = T g and slightly increases when the bath 
temperature decreases below T g . One would instead expect a fictive temperature to 
remain stuck to T g in a (mean-field) model in which the glass transition is sharp. 
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There are other mean-field models such as the Sherrington-Kirkpatrick spin-glass 
model in which there are infinitely many different effective temperatures. The lowest 
aging-regime temperature appears discontinuously (T (u,t w ) > T g ) as one crosses T g 
and can be shown || to decrease with decreasing temperature. Another example of 
this kind is the model of a particle moving in a random potential with long-range 
correlations [f2~4] . p£}l . 



In all cases we know of domain growth models |23| one obtains X(C) — > for t w — > oo, 
in the whole low temperature phase. Hence To(uj,t w ) — > oo in the aging (coars- 
ening) regime for all heat bath temperatures below the ordering transition. This be- 
haviour also holds for certain extremely simple disordered systems such as the spherical 
Sherrington-Kirkpatrick spin-glass lj46j , the toy domain growth model we used in this 



paper. 

It is important to remark that these simple examples have the (sometimes unrealistic) 
feature that nothing depends permanently upon the cooling procedure. One expects, how- 
ever, that more refined models that go beyond the mean-field approximation will capture a 
cooling rate dependence that will also become manifest in the effective temperature. 

Another attempt to identify and relate a microscopic effective temperature to the Active 
temperature of glasses was put forward by Baschnagel, Binder and Wittmann [K71 in the 
context of a lattice model for polymer melts. They have pointed out that in this model 
the usual FDT relation between the specific heat with the energy fluctuations is broken at 
low temperatures and have tried to identify the FDT violation factor with an expression 



they propose for an internal temperature. The internal temperature defined in Ref. f47 
has, however, the unpleasant property of not reducing to the bath temperature in the high 
temperature phase. 

In Ref. [|IH] Franz and Ritort have also discussed the possibility of relating the FDT ratio 
with an effective temperature, in the particular case of the Backgammonn model. They have 
compared the the value of X(t, t w ) for finite times t, t w with the temperature arising from an 
adiabatic approximation they used to solve the model. The result is negative, in the sense 
that T/X(t,t w ) does not coincide with the "adiabatic temperature". 

These examples suggest that the phenomenological Active temperatures act essentially 
as parameters for describing an out-of-equilibrium "equation of state" while the effective 
temperature we have discussed plays a role closer to the thermodynamical one. 



VII. DISCUSSION, EXPERIMENTAL PERSPECTIVES AND CONCLUSIONS 

Although our discussion has been biased by the models we can solve at present, we 
feel that the concept of the effective temperature that we discussed should be relevant for 
many systems with small energy flows. Indeed, the key observation we make is that this 
temperature actually controls thermalization and heat flows within a time-scale. Therefore, 
the effective temperature can be a starting point for the investigation of the thermology and, 
hopefully, the thermodynamics of out-of-equilibrium systems with small energy flows. 

In order to use this idea as a guiding concept for the planning of experiments, one has 
to take some care: obviously, waiting times can be large but not infinite in experiments, 
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and stirring rates can be small but not infinitesimal. Moreover, the models that we have 
explicitly discussed do not exhibit dependence on the cooling history of the sample. Now, 
one would expect in general the effective temperature to depend on t, t w and, e.g., on the 
cooling rate. One can hope to catch this aspect in more refined models. 

One should also pay attention to the equilibration times. The effective temperature 
measured by a small oscillator is related to the FDT violation factor provided that the 
equilibration time t c (Eq.( |3.f OP ) remains much smaller than t w . It sometimes happens that 
an observable O which should thermalize at an effective temperature To{w,t w ) exhibits so 
small fluctuations at this frequency that this time becomes unbearably long. This is the case, 
for instance, of high k modes in the aging time scale for a manifold in a random potential. 
Therefore, although all k modes nominally thermalize at the same effective temperature, 
only low enough k modes can be effectively used to measure it. How low k must be will 
depend on the time scale one is looking at. 

Several numerical and real experiments in structural glasses can be envisaged. For ex- 
ample, one can compare density-fluctuations and compressibility in different length scales 
in order to check if they are equilibrated within a time-scale. Since the low-temperature 
extension of the MCT makes definite predictions on the value of the first non-trivial effec- 
tive temperature appearing as one crosses the transition (see Section VI), this provides a 
concrete ground for a experimental testing of MCT. 

We close by considering the following (slightly Gedanken) spin-glass experiment of Fig. 



10] : Currents are induced in the coil by the magnetization noise of the spin-glass, which is 
in contact with a heat bath. Apart from the interaction with the sample, the L-C circuit 
of coil and capacitor is without losses. This is exactly a realization of the oscillator as a 
thermometer of Section II. From what we know from the time-scales of real spin-glasses 
ZD| , |2"2]] , if the time after the quench is of the order of fO minutes, and the period of the L-C 



circuit is of the order of the second, we are probing (at least partially) the aging regime: the 
temperature (defined as the average energy of the capacitor) should be different from the 
bath temperature. We believe that it would be interesting to return to the magnetization 
noise experiments [21] with the purpose of measuring the effective temperature: this would 
give us, for instance, useful insights into the nature of the spin-glass transition. 
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FIG. 10. An effective temperature measurement for a magnetic system. The coil is wound 
around the sample, which is in contact with the bath. Coil and capacitor have zero resistance 



APPENDIX A: THE TEST MODEL 

Throughout this paper we have used as a test model a spherical disordered model with 
p-spin interactions. The model has been introduced in its purely relaxational version by 
Crisanti and Sommers [4j| as a simple spin-glass model with several advantages, in partic- 



ular, that exact dynamical equations can be written for it in the thermodynamic limit. As 
shown by Franz and Hertz [(5J, these dynamical equations are also those obtained from the 
MCA to the Amit-Roginsky model |J£|. When considered in full generality, the two-time 
dynamical equations correspond to the low temperature extension ||40| , |33| of the simplest 
mode coupling theory for the supercooled liquid phase proposed by Leutheusser [[Kj and 
Bengtzelius, Gotze and Sjolander |]5T| , |34| . It has been also recently shown by Chandra et al. 
that this model is related to a mean-field approach to Josephson junction arrays [52] . 
We consider a system of N variables s = (si, . . . , s^), subject to forces Ff given by 

•'^ E , (Al) 

{jlvjp-l} 



s 



where the couplings are random Gaussian variables. For different sets of indices 
{i,ji, . . . , jp-i} the J's are uncorrelated, while for permutations of the same set of indices 
they are correlated so that 



F/(s>)F/(s) = SiMq) + a^M/N, (A2) 

where q — (s ■ s')/N. In the purely relaxational case, one has f2{q) = f[{q)- We take here 
/ 2 (g) = (1 — D)f[(q), where fi(q) = pq p ~ 1 /2. The couplings jf 1 '"^ -1 are symmetric under 
the permutation i <-> in the purely relaxational case (D = 0). On the other hand, if 
j3j2:.j P -i anc j jm-jp-i are uncorre i a ted, one has D — 1. 
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The dynamics is of the Langevin type: 

di = -F/(s) - fi(t)si - hi(t) + rji(t), (A3) 

where rj is a white noise of variance 2T, fi(t) is a Lagrange multiplier enforcing the spherical 
constraint E^Li s 1 — 1> an d is an external field (usually set to zero), needed to define 
the response functions. 

In Figs. [I] and [| we plot the x{C) curves for the asymmterical and symmetrical p = 3 
model, respectively. In Figs. ^ and |] we plot the auto-correlation decays for the symmetrical 
and asymmetrical versions respectively. Figures [| and ^ show the x{C) curves for the p = 2 
version that is analogous to the O(N) model in D = 3 when iV — > oo. 

In order to check thermalization in this particular example we consider the evolution of 
two such systems with spins s and a, and uncorrelated realizations of disorder J and J' and 
thermal noise r\ and rf, respectively. They are coupled via the term proportional to fi\2 and 
H21 in the Langevin equations 

ds ' 

"T7~ = -Fz J (s) ~ A*ii (*) «i - A*i20i + hi + r]i , (A4a) 
at 

= -KFf{a) - ^22 (t) ^ - /iaiSi + ^ + t?< . (A4b) 

We set = /i 2 i, so that the coupling does not contribute to the stirring. The coefficients 
fin are the Lagrange multipliers for each system. The factor K is introduced to break the 
symmetry between the subsystems. The correlations 



C n (Mw) = £££i(si(t)st(tw)> , C 22 (t,t w ) = jfE? = M(t)<T i (t v )) , 



(A5) 



and responses 

flu(t.O = jjJ2l 1 5(s t (t))/5h s l (t w ) , i? 22 (Mw) = j f l£iS(<Ti{t))/6hZ(t w ) , r . . 
i?i 2 (Mw) = ^Eili5(si(t))Mf(tw) , #21 (Mw) = ^Eili5(si(t))Mf(tw) , 1 J 

precisely satisfy Eqs. ( |5.1a| )-( |5~4D with F given by (|5.21| ). 

The results for the effective temperatures obtained from the numerical integration of the 
exact evolution equations of these systems are shown in Figs. ^ and ||. See main text for the 
discussion. 



APPENDIX B: ENERGY OF THE OSCILLATOR 

In this Appendix we solve Eq. (|3.8| ), we compute ^uj 2 {x 2 (t)), the average potential energy 
of the oscillator. We thus prove Eq. ( |3.14| ) and its form Eq. ( ^.13| ) valid for the stationary 



case. 

Let us define 

rt 



X(w,t) exptiut) = I dt' R(t,t')exp(iut') , (Bla) 
Jo 

C(u,t) exp(ioot) = ( dt' C(t,t')exp(iut') . (Bib) 
Jo 
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If uo^ 1 <C t, we can assume that x(u;,t) and C(u,t) are functions that vary slowly with t, 
thus defining a Fourier component that is "local" in time t. 
In general 

( x 2 (t)) = a 2 C dt' C dt" G(t, t') G(t, t") C(t', t") , (B2) 

Jo JO 

where G(t, t') is the Green function for the oscillator plus the term representing the response 
of the system and C(t,t') is the system's auto-correlation function. Using the definition in 
Eq.( pia|) one can show that the damped oscillator's Green function reads 



*) = 2~IT 2~~ 2 1T( > ( B3 ) 

G(t, = exp f-^y) M^o(t - 0) 0(t - t') . (B4) 

We have here replaced x(t,w) by x"(t, u ) using the fact that a 2 N <C 1. The characteristic 
time t c (t) of the damped oscillator is given by 

= ^i) ' (B5) 

We can now study Eq.(|B"2h by using the above expressions for G(t,t'). After a simple 
change of variables, using causality and the fact that G(t,t') decays exponentially as a 
function of time-differences, Eq.( p2|) can be rewritten as 

/oo rji.) roo rco 
— dr dr' G{t, t - r) G(t, t - r') C(u, t - r) exp(^(r - r')) . 
-oo 27T J —oo J —oo 

(B6) 

The fast exponential decay of the Green function allows us to replace C(u,t — r) by C(u, t). 
Thus, 

(x 2 (t)) = a 2 — G(u, t) G{-cu, t) C(u, t) = a 2 — — 



x 2tt v y v ' v ' ' 7oo 2vr x"(w,t)-x"(-w,*) 
1 1 



^ - cj 2 - a 2 iVx"(-^, t) lu 2 - lj 2 - a 2 N X "(uj, t) 



(B7) 



This integral can be calculated by the method of residues. We can close the circuit on 
the upper complex half plane. Since x"(u;,i) is analytic in the upper half plane, the only 
singularities are the zeroes of the denominators. Assuming that a 2 x" is small, they lie in the 
vicinity of u — ±u> . In fact one can check that only two poles penetrate inside the circuit 
of integration. We obtain therefore 



This is the general result for the temperature. The particular result ( p.!3|) that holds for 
the stationary case is recovered from Eq. (|B8|) by letting C(t, uS) and be independent 

oft. Thus, flBHp reduces to (gTjj). 
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APPENDIX C: SMALL BUT MACROSCOPIC THERMOMETERS 



In this Appendix we show that the thermometric considerations made in Section III do 
not crucially depend on the choice of an oscillator as a thermometer. We use here a small 
but macroscopic thermometer defined by the variables y^ i — 1, . . . , n, and we couple it only 
to the observable O(s) of the system through a degree of freedom x(y). 

Our measurement procedure is as follows: we first thermalize the thermometer with an 
auxiliary bath at temperature T*. We then disconnect it from the bath and we connect it 
to the system through O. If there is no flow of energy between thermometer and system, 
then we conclude that the measured temperature is T*. 

The energy of the thermometer plus its coupling with the system is 

H = E(y) - aO(s)x(y) . (CI) 

The net power gain of the thermometer is then Q(t), given by 

Q = a(Ox) = a d t ,(x(t)0(t'))\ t ,^ t - . (C2) 

We look for the condition that ensures stationarity for the thermometer. The thermometer is 
characterized by a temperature-depedent correlation C x (t, t') = (x(t)x(t')) and its associated 
response R x (t,t'). Using linear response one has 

Oft) = O h (t) + a f alt' Ro(t, t')x(t') , (C3a) 
Jo 

x(t) = x b (t) +a [ dt' RM,t')0(t') . (C3b) 
Jo 

To leading order in a, (Ox) is given by 

(Ox) = af dt" (R x (t, t") C (t", + Ro{t', t") C x (t, t")) . (C4) 
Jo 

Assuming now that T* is such that the thermometer can be considered to be almost in 
equilibrium (note that the coupling a is small and that we have chosen a thermometer that 
is not itself a glass!), we obtain 

Q = a 2 J*dt> R x (t-t') { ^ C °^ t ' ) -T*R (t,t')\ . (C5) 

The condition for having no flow is then that the average of the parenthesis in the integral 
is zero. The weight function for this average is R x (t — t') which contains the characteristic 
time of the thermometer. 



APPENDIX D: DEFINITION OF TIME-CORRELATION SCALES 

In this appendix we review briefly the definition of correlation "time-scale" introduced 
in Ref. H] for a correlation function that depends nontrivially upon two times. 
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Given a correlation function C(t,t w ), which we assume normalizable in the large-times 
limit C(t,t) — > Coo > 0, we consider three increasing times t\ < < t 3 , and the limit in 
which they all go to infinity, but in a way to keep C{t^,t2) = a and C(t2,ti) = b constant. 
We thus define the limit 

lip C(* 3 ,ii) = /M). (Dl) 

C(t 3 ,t 2 ) = a , C(t 2 ,tl)=i> 

The mere existence of the limit "triangle relation" / has extremely strong consequences: 
considering four times one can easily show that / is associative 

f(a,f(b,c))=f(f(a,b),c). (D2) 

The form of an associative function on the reals is very restricted, and a classification of all 
possible forms can be made ||. 

It is sometimes convenient to work with the "inverse"/ of / defined as 

f(a,b)=c f{a,c) = b. (D3) 

We can now define a correlation scale in the following way: given two values of the 
correlation at large times C(t, t w ) = C 1 and C(t,t' w ) = C 2 , t > t' w > t w and C 1 < C 2 , they 
are in a different correlation scale if 

f(C 2 ,C 1 )=C\ (D4) 

and are in the same scale otherwise. In other words, the time it takes the system to achieve 
C 2 is negligible with respect to the time it takes to achieve C l . 



APPENDIX E: SOLUTION OF THE TWO COUPLED-MODES EQUATION 

Using the separation (|5.71) - ([5.14|) and the mode-coupling approximation ( p73|) -( pl4|) we 
obtain a similar separation for D a b and E Q ^: 

Dab DT (t - t w ) = F ab {C FDT {t - t w )) , (Ela) 
Kb DT {t -*w) = £ F abyCd {C FDT {t - t w ))R F cd DT {t - t w ) , (Elb) 

c.d 

where C FDT stands for the set C F b DT (t — t w ). One similarly obtains, in the aging regime, 

D ab (t,t w ) = F ab (C(t,t w )) , (E2a) 

£af>(Mw) = / J Fab,cd(C(t, tyr))Rcd(t, tw) • (E2b) 

We can now write two coupled sets of equations, valid in the quasiequilibrium regime and 
in the aging regime, respectively. For t — t w finite, and large times t > t w , we have: 
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at 



E 



FDT 



oo I ^ ac 
H'ac ' 7 



(0) 



T 



Y,D E ac A Cf b A + M t 



DO 

ab 



4?" 

where /i~ = linn^ /i ac (t) and 



dt" D F ac DT (t-t")- cb 1 v 



0*. 



In the aging regime, for t > t w we have 



dC a b(t, t v 



dT~ 



E 



A*ac(*) + 



^(0) - D. 
T 



EA 
ac 



Ccb{t, t w ) + D ac (t, t v 



(E3) 



(E4) 



^ DT (Q) - cg A 

T 



+ E / W d *" ^«c(t, O ^*(*w, *") + E / dt " Vac(t, t") C cb (t", t w ) , (E5a) 



E + ^ — — 

c£ OT (o) - ci A 



T 



Rcb(t,t w ) + E / ^ac(^,^") Rcb(t",t w ) 

(E5b) 



Equation (|E4|), for given is very similar to the high-temperature Mode-Coupling equa- 
tions ||34|| , and can be solved in the same way. An asymptotic solution for the aging regime 
can be obtained by using the generalization to more than one mode of the ansatz in Ref. || 
( |5.15| ). The derivative terms in Eqs. ( |E5| ) can be then dropped provided that X n ^ and 
X 2 2 7^ 0, a fact to be verified a posteriori. We shall find in this way two different solutions 
for Eqs. (|E5|) . 

In the unthermalized case Xn and X22 are different from zero, and possibly different from 
each other, while X 12 = X 2 \ = 0. It is then easy to see that Eqs. (|E5|) become effectively 
uncoupled in this regime, and can be solved || as two separated one-mode equations, with 
the ansatz 



C aa {t,t w ) — C aa (h aa (t w ) / h aa (t)) , 

Raa{t,t w ) = — — — (t,t w ) , 
Cl2[t,t w ) = C 21 (t,t w ) = . 

In the thermalized case, one assumes 

Xn = X22 — X\2 = X21 = X 7^ 
hu = h 22 = hi2 = h 2 i = h . 



Making the change of variables 



h(t w ) 
hit) 



h(t") 



(E6a) 
(E6b) 
(E6c) 



(E7a) 
(E7b) 



(E8) 
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one finds that only the dependence on A survives in the equations, and that they reduce 
to a set of four (instead of eight) consistent equations for the correlations (see [|,^] for a 
systematic approach). 

In this way, both for the thermalized and the unthermalized case one can obtain X a b, 
M a b and C^ b A . One then has to check that Xn and X22 7^ 0. If this is not the case the 
equations become identities, and one cannot anymore neglect the derivative term. One has 
therefore to use a more refined long time limit. These values have to be substituted in (|E3|), 
in order to complete the solution in both regimes. 

Let us remark here that the problem of selecting the functions h ab remains open. This is 
an asymptotic matching problem in a non-local equation, and does not appear to be easily 
solvable. 

We have thus found the long-time limit of the correlations and responses. If there 
is more than one asymptotic solution (even modulo h), we don't know for the time being 
which asymptotic form is selected by the unique solution of the evolution equations, without 
resorting to explicit numerical integration. 

In short we have: 

1. Unthermalized aging regime. 

X 12 ->0, X 21 ->0, X U ^X 22 . (E9) 

We have therefore 

flu(t,*w) = ^^, %(Mw)4f ' ( E1 °) 



and 



£n(«.«.> = ^. *W«.0 = ^. (Ell) 

2. Thermalized aging regime. 

X n = X12 = X 2 i= X 22 = X , (E12) 
O = f Va,6. (E13) 

where the aging time-scales "lock in", i.e., there is the same function h(t) for all a, b, 
such that 

C ab (t,t w ) = C ab (h(t w )/h(t)) . (E14) 

This property can also be stated by saying that as t,t w — >• 00 a plot of Cn(t,t w ) vs. 
C22{t,t w ) yields a single smooth curve: i.e., that there is a function Ti(C) such that 

Cn(t,t w )=H(C22)(t,U) . (E15) 
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